Non-Abelian operations on Majorana fermions via single charge control 
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We demonstrate that non-Abelian rotations within the degenerate groundstate manifold of a set 
of Majorana fermions can be realized by the addition or removal of single electrons, and propose 
an implementation using Coulomb blockaded quantum dots. The exchange of electrons generates 
rotations similar to braiding, though not in real space. Unlike braiding operations, rotations by 
a continuum of angles are possible, while still being partially robust against perturbations. The 
quantum dots can also be used for readout of the state of the Majorana system via a charge 
measurement. 



Elementary excitations of topological materials can 
have unusual properties, such as statistics different from 
that of fermions or bosons. The most interesting possibil- 
ity is non-Abelian statistics which is believed to realized, 
for example, in the 5/2 fractional quantum Hall system 
and in topological superconductors, where the quasipar- 
ticles are Majorana fermions [TJ. Following the original 
suggestion by Kitaev[2], recent theoretical achievements 
have shown that topological superconductors should be 
realizable in semiconductor or metallic systems with the 
right combination of spin-orbit coupling, induced super- 
conductivity and applied magnetic field [3HT3]. which has 
brought new optimism into the search. Several methods 
to observe the non-Abelian nature of particle exchange 
have been suggested. For fractional quantum systems, in- 
terferometry of paths (TUE3 an d Coulomb blockade peak 
spacings|16j sensitive to the number of enclosed quasipar- 
ticles have been analyzed, and the former recently inves- 
tigated experimentally [1 7 . The non-local nature com- 
bined with Coulomb interactions has also been exploited 
theoretically [TB]. For semiconductor wire systems, ex- 
change of particles in one-dimensional network structures 
was recently proposed as a direct way to observe non- 
Abelian features of Majorana bound states [T5]. 

Moreover, non-Abelian quasiparticles have been sug- 
gested as a basis for topological quantum computing [20]. 
where a computational step is done by physically ex- 
changing the positions (braiding) of the quasiparticles, 
thus performing a unitary rotation in the degenerate 
groundstate manifold. In principle, the exchange opera- 
tion only depends on the topology of the exchange and 
the system is therefore argued to be robust against lo- 
cal perturbations [21 [T3] . However, braiding of Majorana 
fermions has to be supplemented with gates that are not 
topologically protected in order to give a set of universal 
operations. Several methods have been put forward, in- 
cluding bringing the quasiparticles together for a certain 
amount of time [ST], and very recent proposals combine 
topological qubits and conventional qubits [2"2Tf2~i] for im- 
plementation of quantum gates. 

In this Letter, a different method for performing opera- 



tions on a set of Majorana bound states is analyzed. The 
individual operations are adiabatic tunnel processes of a 
single electron from a Coulomb blockaded quantum dot 
coupled to one or two Majorana modes. The processes 
result in rotations of the groundstate manifold, reminis- 
cent of actual physical exchange of the particles. Like 
real-space braidings, the tunnel-braid operations are ro- 
bust against dephasing. However, they are sensitive to 
electrical noise on the tunnel amplitudes, and in this re- 
spect not topologically protected. The electrical control, 
on the other hand, allows rotations with tunable angle. 
Moreover, to achieve the maximal protection of the tun- 
nel braids, a certain phase difference (controllable by a 
flux) between the two Majorana bound states is needed. 

The system we have in mind is a network of semi- 
conducting wires or a two-dimensional electron gas with 
spin-orbit coupling, induced superconductivity and ap- 
plied magnetic field, and in addition a number of quan- 
tum dots, see the sketch in Fig. [2]Ja). For parameters 
where the superconductor is in the topological phase, Ma- 
jorana bound states (MBS) exist at the ends of the wires 
[7HT0] . They are zero energy solutions to the Boguliubov- 
de Gennes equations and have the general form 

7* - / dr(/i* t + #*? + 3**4. + ■ W 

The MBS operators obey ji — 7J and jf — 1. Here 
fi and gi are functions of the spatial coordinate r and 
^cr = 4 r (T (r) is the electron field operator. The Majo- 
rana bound states are tunnel coupled to quantum dots 
in Coulomb blockade regime, i.e. the charge state on a 
dot is restricted to be either N ox N +1 electrons, with N 
some arbitrary number. The level spacing of a dot is sup- 
posed to be much larger than temperature, so that only 
one quantum state is involved. Furthermore, because a 
large magnetic field is applied to induce the topological 
superconducting phase, it is assumed that spin degener- 
acy of the dots is broken, so that we only have to consider 
a single spin direction, say, spin up. The dot potential 
energies can be controlled by electrostatic gates. 

Projection onto Majorana bound states {7, } coupled to 
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even and odd parity the Hamiltonian matrix is 
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FIG. 1. (Color online) A single Majorana state (Ml) cou- 
pled to a quantum dot (D) by a tunnel coupling Vi. The Ma- 
jorana state Ml is combined with another Majorana bound 
state, M2, to form a fermion, M12. Starting in the ground- 
state, the operation Pi takes the energy e (using the gate (G)) 
from negative to positive values causing an adiabatic transi- 
tion from a full dot to an empty dot. The electron is added 
to the superconductor, thus inverting the parity of the Ma- 
jorana system. Because the total even (red) and odd (blue) 
parity states are degenerate the inversion is independent on 
the details of the transition. 



an electron level (described by a spin-up electron annihi- 
lation operator c) leads to an effective low-energy tunnel 
Hamiltonian 
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where Vi is the Hamiltonian overlap between the electron 
wavefunction, cj>, and the spin-up electron component of 
the MBS Vi = (fi\H\cf>). This matrix element will be 
important for the actual rotation done to the Majorana 
system and should therefore be controllable. 

First, we consider the situation with a single MBS 71 
coupled to a single dot and study the transition of the 
groundstate as the occupancy of the dot is changed adi- 
abatically. Two Majorana fermions are needed to define 
a usual (Dirac) fermion. Using Majorana states 71 and 
72 as basis, a fermion M12 with annihilation operator 
d = (71 + 272) /2 is defined. The two eigenstates of <ffd 
are denoted |0)mi2 and 1 1) A/12 ■ The Hamiltonian for a 
dot level c connected to a Majorana mode 71 = d + eft is 
then 



Hi = £c f c+ (v*^ - vic)(d + dft), 



(3) 



where e is the dot level energy, measured relative to the 
chemical potential of the superconductor. The Hilbert 
space contains 4 states, and the Hamiltonian is block 
diagonal by conservation of the total parity. For both 
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with the basis for even and odd total parity cases being 
{|0) d |0)mi2,|1)d|1)mi 2 } and {|0) d |1)mi2, |l) D |0) M ia}, 
where 0(1) represent an empty(full) fermion state, and 
|-)d denotes the dot state. 

Now an operation Pi that adiabatically changes e/v\ 
from —00 to +00 is defined, see Fig. [I] If the original 
state of the Majorana system is \i)m = c*|0)mi2+/3|1)mi2 
and the dot state is |1)d, the state of the dot plus Ma- 
jorana system is at any gate potential e given by the 
superposition (up to an overall dynamical phase factor) 
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where ui&(e) = Ea(e) and E = e/2 — ^/{e/2) 2 + v\. The 
dot and the Majorana system is thus entangled during 
the operation, but not at the end of the operation where 
e/vi — > 00 and b — > 1. Therefore, the operation P\ gen- 
erates an inversion of the occupation of the M12 fermion, 
which in basis independent notation|25] is expressed as 



Pi 



IVm !-> 7iIVm- 
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If a number of Majorana bound state are connected 
to quantum dots, it is thus possible to manipulate the 
Majorana system by repeated applications of P opera- 
tions, yielding a new state 71 • • . 7 to K)m- However, be- 
cause {7i,7j} = 2Sij this only provides a finite number 
of operations. Interestingly, two consecutive operations 
P\Pi corresponds to a real-space process where one Ma- 
jorana bound state is rotated around the other [57]. 

A natural question is how sensitive the operations Pi 
are to decoherence of the dot charge. To investigate this, 
we add environmental degrees of freedom that couple to 
the dot charge. At beginning of the P\ process the ini- 
tial state is = |l)D|i)M|?i)cnv, where | 7i) cnv is the 
initial environment state. When expanding in powers of 
the tunneling Hamiltonian, the time-evolution operation 
can be collected in powers of 7! as U (T, 0) = U a + 7if7i. 
If the transition is done slowly enough to guarantee that 
the occupancy of the dot changes by one (can be checked 
by measuring the charge on the dot) only the part jiUi 
survives, and Eq. ^ holds, regardless of the dot's cou- 
pling to the environment. 

As mentioned, the single MBS processes in Eq. ^ only 
allow a limited set of operations. A richer set is possible 
if the dots are connected to two Majorana modes, as in 
Fig. [5Ja) . Using the same basis as above, the even and 
odd sector Hamiltonians for MBS 71 and 72 coupled to 
a dot (Dl in Fig. 2) are in this case 
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FIG. 2. (Color online) (a) An one-dimensional array of Majo- 
rana states (Ml,...,Mn) coupled to quantum dots (D1,D2,...) 
in the Coulomb blockade regime. Each dot is tunnel coupled 
to two Majorana states with tunnel barriers (controlled by 
the gates adjacent to the plunger gates (G1,G2, ...)). Chang- 
ing the occupancy of a dots by one electron creates the unitary 
rotations Pij. (b) The groundstate energy of one dot coupled 
to two Majorana bound states, with |«i| = |ti2| for even (red) 
and odd (blue) total parity of a dot and its two connecting 
MBS. Even and odd cases are degenerated for tpi — 2nir, 
which makes the operations partially protected. The full 
and dashed lines are for e/|«i| = and 2, respectively. 



where v CV cn/odd — v i T * w 2- In general, the matrix el- 
ements Wcvon and u dd are different and the degeneracy 
between the even and the odd cases is lifted. The ground- 
state energies are (see Fig. 2(b)) 
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where v 2 = |«i| 2 + \v 2 \ 2 and (pi = 2Arg(u 1 /u 2 ). The 
phase different is controlled by the flux $i and up to 
constant phase shift, we can write ipi — $i/$0i with $i 
being the flux in loop 1 and $o = h/2e. Note that Eq. Q 
is 47T periodic in ip-\ }26j. 

Again, consider an adiabatic process that transfers an 
electron from a dot but this time to the two Majorana 
states 71 and 72. Unlike the case with a single MBS, 
the resulting rotation of the Majorana system is in gen- 
eral not independent on the time spend in the adiabatic 
process, because of the energy difference in Eq. see 
Fig. gb). The de generacy is restored only when the 
phase difference is ip± — 2nit (n integer) , which therefore 
requires tuning the magnetic flux $1 (Fig. [2^a)). At this 
degeneracy point, V1/V2 is real which allows the Hamil- 
tonian to be written as 



where a new Majorana operator is defined 
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(10) 



and where a common phase is absorbed into the dot 
electron operator c = cexp(zArg(vi)). Thus, since the 
Hamiltonian ([9| has the same form as ([3|, a dot coupled 
to two MBS reduces (at the degeneracy point) to a dot 
coupled to a single Majorana state 712. The conclusion 
from above therefore also carries over: by adiabatically 
changing the electron number of the dot, the following 
rotation is performed 
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To understand the rotations that can be generated 
by repeated applications of Pi 2 (with different ratios 
I v i l v i I ) i we use the following Pauli matrixes acting on the 
two level system spanned by 71 and 72 : a x — 71 , a y — 72 , 
and o z — — «7i72- In this language, the operation P12 
makes a 7r-rotation around an axis in the x-y plane, but 
other rotation angles around lines in the x-y plane can- 
not be done. In contrast, when applying a pair P12P12 = 
(u r fi+v~f2){u'~fi+v'j2) — (uu' +vv') + i(uv' — vu')a z aro- 
tation around the z-axis with tunable angle is performed. 
A braid operation also rotates around the z-axis, but by 
an angle restricted to 7r/2. Instead, using four MBS and 
the even-parity subspace to define a qubit [2T] . a universal 
set of single qubit rotations is in fact generated by pairs 
of P operators. Again, Pi2-P{ 2 IS a rotation around the 
z-axis (in the basis {(00), (11)} defined below), whereas 
P23P23 now gi yes a rotation around the x-axis with con- 
trollable angle. However, two-qubit gates are not possi- 
ble in the latter representation using P-operators only, 
see Refs. [TJ I21j for more on quantum computing imple- 
mentations. 

A special and illuminating case is when the dots couple 
with equal strength to two MBS (|vi| = \v2 \ in Eq. (10)), 
which results in operators F s ; = ^qili + Ij+i) acting on 
nearest neighbors. They are related to braid operators 
Bl = 75 + 7i+i7i)E3 by B t = F ai = H+iFi. The F, 
operators fulfill Ff — 1 and 



FiFj = —FjFi, \i-j\>l 
FiFi + \Fi = —Fi+iFiFi+i, 



(12a) 
(12b) 
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which differs by a minus sign from the relations defining 
the braid group|l]. As a side remark, Fi form a projective 
representation of the permutation group [28 . 

To demonstrate the non-Abelian nature of the tunnel- 
braid operations, consider now an explicit example with 
four Majorana states and three dots. The state of the 
superconductor is initialized by tuning the dots and the 
magnetic field to fuse Majorana pairs (1,2) and (3,4) 
and letting them relax. The initial state is 1 00) = 
|0)mi2|0)m34, referring to the occupation of the fermions, 
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d\ = (7i+«72)/2 and d 2 = (73 + 174)72. We will consider 
applications of pairs of Fi and hence restrict to the sub- 
space of even parity, spanned by 1 00) and |11) = d\d\ |00). 
The possible unitary transformation are given by 

(FxF 2 ) eveii = [(F 2 F 3 ) cv J T = ^{ \ 7 ) ' ( 13 ) 
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and {FiF 3 ) cvcn — a x (up to phase factors). Other per- 
mutations can be deduced from FiFj — \FjFj\~ 1 . Appli- 
cation of F1F2 or F 2 Fi gives 

F 1 F 2 |00) = -^(|00) + |11)), (14a) 

F 2 Fi\00) = — (|00}+i|ll». (14b) 
v 2 

For one sequence the resulting state is an eigenstate of 
a x and for the other it is an eigenstate of a y . How- 
ever, expectation values of neither of these operators are 
easily measured using quantum dots, capable of measur- 
ing charge only. Initial application of F 2 F% and F 3 F 2 in 



an occupation measurement and this yields 

F 1 F 2 F 2 F 3 \Q0) = |11>, 
FaFiF 3 F 2 \00) = |00). 



(15a) 
(15b) 



The last sequence can be implemented in the follow- 
ing way: the potentials on the dots Dl, D2, and D3 in 
Fig. [2]j a) are successively increased from positive to neg- 
ative voltages, thus emptying the dots, and then D2 is 
filled again by tuning the D2 potential back to positive 



voltages. The sequence (15a) is done in the same way, 



except in a different order. 

Finally, we discuss how the state of the Majorana sys- 
tem can be read out. Several methods have been pro- 
posed to measure the state of coupled Majorana modes, 
including observing changes in the current-phase relation 
of a Josephson junction with two Majorana states in the 
loop [2 [53], or by interferometry[TJ [15]. In the present 
setup it is natural to use the quantum dots, which is 
indeed possible. By choosing the phase difference be- 
tween two MBS, so that the fused even or odd states are 
split in energy, the occupancy can be read from the adi- 
abatic curves similar to the one in Fig. [l] because the 
off-diagonal elements in Eq. Q, unlike in Eq. Q, differ 
for the even and odd cases. Maximal visibility is achieved 
for ipi = 7r, where the two off-diagonal matrix elements 
are \v±\ ± \v 2 \ for even/odd, respectively. 

In conclusion, a scheme for manipulation of the state 
of a set Majorana fermions has been proposed. It allows 
for demonstration of the non-Abelian nature of the quasi- 
particles, albeit not by real-space exchanges. Instead the 
exchanges are performed via removal or addition of real 
electrons, made possible by Coulomb blockade. 
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